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Preface 


This  thesis  presents  the  results  of  an  investigation  into  the  feasibility  of  using 
frame  selection  and  bispectrum  phase  reconstruction  as  a  post-  processing  technique 
for  Air  Force  speckle  imaging  systems.  The  key  result  is  a  frame  selection  and  bis¬ 
pectrum  phase  reconstruction  algorithm  that  can  be  immediately  implemented  at 
Air  Force  space  surveillance  sites  to  improve  their  mission  effectiveness.  I  hope  that 
this  research  effort  will  contribute  to  our  nation’s  defense  by  enhancing  the  effective¬ 
ness  of  the  Air  Force’s  space  surveillance  capabilities.  Furthermore,  I  hope  that  this 
thesis  will  inspire  future  AFIT  students  to  pursue  research  in  this  challenging  and 
vital  area. 


Elizabeth  A.  Harpold 
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Abstract 

Frame  selection  using  quality  sharpness  metrics  have  been  shown  in  previous 
AFIT  theses  [25]  [5],  to  be  effective  in  improving  the  final  product  of  images  obtained 
using  adaptive  optics.  This  thesis  extends  this  idea  to  uncompensated  speckle  image 
data.  Speckle  image  reconstruction  is  simulated  with  and  without  frame  selection. 
Speckle  images  require  the  processing  of  hundreds  of  data  frames.  Frame  selection  is 
a  method  of  reducing  the  amount  of  data  required  to  reconstruct  the  image  without 
significant  loss  in  quality.  A  collection  of  short  exposure  image  data  frames  of  a  single 
object  are  obtained  and  sorted  based  on  sharpness  metrics.  Only  the  highest  quality 
frames  are  retained  and  processed  to  form  the  final  image.  The  fraction  of  the  frames 
that  are  retained  for  processing  is  user  defined.  The  phase  spectrum  is  reconstructed 
using  the  bispectrum  technique.  The  benefits  of  frame  selection  for  point  (star) 
sources  and  extended  (satellite)  sources  are  examined  by  comparing  composite  image 
data  with  and  without  frame  selection.  The  resulting  power  spectrum  is  evaluated 
through  the  SNR  gain  measurements,  and  the  resulting  phase  spectrum  is  evaluated 
by  measuring  the  phase  error  between  the  data  from  the  composite  image  and  the 
data  from  the  true  object.  In  both  cases,  the  results  show  that  frame  selection  does 
not  improve  the  power  or  the  phase  spectrums.  For  point  sources,  results  show 
frame  selection  causes  slight  decrease  in  performance.  For  extended  sources,  the 
change  in  performance  is  insignificant.  However,  frame  selection  does  offer  a  means 
for  data  reduction  without  significantly  reducing  performance  in  a  wide  variety  of 
target  brightness  levels  and  atmospheric  turbulence  conditions. 
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THE  ROLE  OF  FRAME  SELECTION  AND  BISPECTRUM 


PHASE  RECONSTRUCTION  FOR  SPECKLE  IMAGING 
THROUGH  ATMOSPHERIC  TURBULENCE 

L  Introduction 

1.1  Motivation 

The  Air  Force  has  a  mission  to  monitor  the  status  and  orbit  of  all  space  objects. 
One  means  of  accomplishing  their  mission  is  through  the  use  of  high  resolution 
optical  imagery  obtained  from  ground-based  observatories.  All  observations  from 
ground  based  observatories  are  corrupted  by  atmospheric  turbulence  effects.  The  Air 
Force  employs  various  methods  for  compensating  for  turbulence  effects  on  imaging 
systems.  Adaptive  optics  uses  a  wavefront  sensor  and  a  deformable  mirror  to  help 
compensate  for  some  of  the  atmospheric  distortion.  Conversely,  speckle  imaging  is 
a  post-processing  technique  which  uses  short  exposure,  uncompensated  images.  It 
is  often  observed  that  for  any  given  set  of  short  exposure  images,  there  are  some 
images  that  are  simply  better  than  others  [4,  5,  7,  18,  22,  23,  25].  Frame  selection  is 
based  on  the  characteristics  of  the  individual  short  exposure  image  and  its  associated 
spectrum.  It  is  a  three-step  process  for  choosing  a  subset  of  images  for  further 
processing. 

1.  Calculate  a  frame  quality  metric  value  for  each  frame  in  the  data  set. 

2.  Sort  the  images  in  the  data  set  from  highest  to  lowest  based  on  the  quality 
metric  value. 

3.  Select  the  highest  value  images  for  processing. 
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Previous  work  at  AFIT  has  shown  that  frame  selection  does  indeed  improve  the 
quality  of  compensated  space  images  [24].  During  the  summers  of  1993  and  1994, 
Captains  Craig  Stoudt  and  Steve  Ford  demonstrated  that  frame-selection  can  im¬ 
prove  the  performance  of  adaptive  optics  imaging  systems  when  turbulence  condi¬ 
tions  limit  the  image  spectrum  SNR,  but  not  when  measurement  noise  limits  the 
image  spectrum  SNR  [5,  25]. 

The  primary  goal  of  this  research  is  to  extend  frame  selection  to  speckle  imag¬ 
ing.  Bispectrum  phase  reconstruction  is  incorporated  into  a  new  atmospheric  simu¬ 
lation  technique  developed  at  AFIT.  The  power  spectrum  SNR  and  phase  spectrum 
error  are  analyzed  to  determine  the  effects  of  various  target  brightness  levels  and  tur¬ 
bulence  conditions  on  speckle  imaging  systems.  The  goal  of  frame  selection  research 
is  to  devise  a  method  for  sifting  through  large  sets  of  image  data  frames  and  selecting 
the  fewest  number  of  the  best  frames  that  maintains  or  improves  the  quality  of  the 
final  image.  An  effective  implementation  of  frame  selection  depends  on  the  size  of 
the  data  base  to  be  processed.  For  extremely  large  data  bases,  a  computer  frame 
selection  algorithm  containing  a  simple  but  effective  quality  metric,  may  be  used  to 
effectively  reduce  the  amount  of  data.  The  selected  subset  may  then  be  processed 
using  a  bispectrum  phase  reconstruction  algorithm. 

1.2  Problem  Statement 

This  thesis  examines  the  role  of  frame  selection  for  uncompensated  speckle 
images,  and  measures  the  performance  of  frame  selection  by  using  bispectrum  phase 
reconstruction  and  phase  error  analysis. 

For  each  case  study,  the  following  user  parameters  are  varied: 

•  Sharpness  metric 

•  Frame  Selection  Rate, 

and  the  following  environmental  parameters  are  varied: 
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•  Vo  (Pried  seeing  parameter  [6]) 

•  (Visual  Magnitude,  or  brightness  levels). 

In  comparing  the  cases,  the  results  are  in  the  form  of 

•  SNR  plots 

•  SNR  Gain  plots 

•  Phase  Error  Plots 

•  Comparative  Phase  Error  Difference  Plots 

•  Deconvolved  Images 

1.3  Approach 

This  thesis  addresses  the  problem  by  investigating  frame  selection  performance 
under  a  variety  of  conditions.  The  power  and  phase  spectrums  of  the  reconstructed 
image  are  measured  to  determine  relative  performance  with  and  without  frame  se¬ 
lection.  To  achieve  these  goals,  the  following  steps  are  undertaken: 

1.  Modified  existing  code  to  account  for  photon  shot  and  CCD  read  noise  sources. 

2.  Modified  code  to  include  bispectrum  phase  reconstruction. 

3.  Developed  a  phase  error  metric  for  comparing  performance. 

4.  Developed  a  matrix  of  computer  runs  to  evaluate  performance  based  on  a  set 
of  key  optical  parameters. 

5.  Developed  a  computer  program  for  analyzing  the  data. 

6.  Performed  runs 

7.  Analyzed  data 

8.  Drawn  conclusions 
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This  research  eifort  relies  on  computer  simulation.  Michael  C.  Roggemann 
developed  a  simulation  package  known  as  HYBISPSIM  which  simulates  imaging 
through  turbulence  and  calculates  statistical  quantities  used  for  speckle  imaging  [25]. 
Another  Roggemann  code,  referred  to  as  HYSIM,  also  simulates  imaging  through 
turbulence,  but  does  not  compute  speckle  imaging  quantities.  In  a  separate  research 
effort.  Captain  Craig  A.  Stoudt  modified  HYSIM  to  incorporate  statistical  frame 
selection  and  a  charge  coupled  device  (CCD)  camera  model  [25].  This  modified  code 
was  renamed  HYSIMS  [25].  This  thesis  effort  combines  HYBISPSIM  and  HYSIMS 
and  modifies  HYSIMS  to  do  speckle  imaging.  A  second  simulation  program,  called 
BIPHASE,  is  created  from  a  series  of  FORTRAN  routines,  also  provided  by  Rogge¬ 
mann,  which  reconstructs  the  phase  from  the  bispectrum,  measures  the  phase  error, 
implements  the  pseudo- Weiner  filter  [21],  and  reconstructs  the  image. 

1.4  Scope 

The  scope  of  this  thesis  is  restricted  to  the  following  topics: 

•  Extend  statistical  frame  selection  work  to  speckle  imaging 

•  Extend  current  AFIT  simulation  capabilities  to  include  the  bispectrum  tech¬ 
nique,  which  is  insensitive  to  tilt  error. 

•  Investigate  alternate  sharpness  metrics  for  frame  selection. 

r 

1.5  Chapter  Outlines 

The  following  is  a  brief  synopsis  of  the  information  found  in  each  chapter  of 
the  thesis. 

1.5.1  Chapter  2.  This  chapter  presents  an  overview  of  the  problems  as¬ 
sociated  with  processing  space  image  data,  including  atmospheric  conditions  and 
equipment  hmitations.  It  also  introduces  the  concept  of  frame  selection  as  a  strat¬ 
egy  for  correcting  for  some  of  these  adverse  affects.  Chapter  2  also  details  speckle 
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image  processing,  including  the  theory  of  speckle  interferometry  and  bispectrum 
phase  reconstruction. 

1.5.2  Chapter  3.  This  chapter  lays  out  the  approach  to  simulating  the 
frame  selection  process  and  the  development  of  bispectrum  experiments. 

1.5.3  Chapter  4-  This  chapter  presents  the  results  of  the  series  of  frame 
selection  case  studies  for  point  sources  and  extended  sources. 

1.5.4  Chapter  5.  This  chapter  states  the  conclusions  based  upon  the 
results  presented  in  chapter  4.  It  also  concludes  with  recommendations  to  the  user 
concerning  the  implementation  of  frame  selection  as  it  applies  to  speckle  image  data 
and  recommendations  for  further  research. 

1.5.5  Appendix  A.  This  appendix  contains  plots  of  data  obtained  in  the 
point  source  experiments. 

1.5.6  Appendix  B.  This  appendix  contains  plots  of  data  obtained  in  the 
extended  source  experiments. 

1.6  Summary  of  Key  Results 

This  investigation  demonstrates  that  frame  selection  for  point  and  extended 
sources,  such  as  satellites,  does  not  significantly  improve  nor  degrade  the  power 
or  phase  spectrums  of  the  reconstructed  image.  This  thesis  demonstrates  that  the 
primary  benefit  of  frame  selection  to  speckle  imaging  is  data  reduction.  Data  can  be 
reduced  by  as  much  as  50%  with  no  loss  in  image  quality. 

1.7  Conclusion 

Atmospheric  turbulence  is  a  major  impediment  to  Air  Force  space  surveillance. 
Frame  selection,  although  incapable  of  fully  compensating  for  these  effects,  offers 
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an  effective  strategy  for  reducing,  at  least  by  50%,  the  data  required  in  the  post 
processing  of  speckle  image  data. 


11.  Background 


2.1  Introduction 

This  chapter  provides  the  background  necessary  to  understand  some  key  el¬ 
ements  of  measuring  and  processing  Air  Force  space  surveillance  and  astronomical 
image  data.  The  first  section  explains  some  of  the  terms  associated  with  speckle 
imaging.  Next,  the  problem  of  atmospheric  turbulence  and  how  it  adversely  affects 
the  quality  of  these  images  is  discussed.  Then,  a  method  for  correcting  for  some 
of  the  turbulence  effects  is  detailed  in  three  parts.  First,  speckle  interferometry,  a 
method  for  obtaining  the  object  power  spectrum  is  explained.  Then  the  bispectrum 
technique  of  phase  reconstruction  is  discussed.  Finally,  frame  selection,  a  method 
of  selecting  and  discarding  data  from  the  detected  image  of  the  object  is  outlined. 
Together,  all  three  elements  provide  an  effective  post-processing  technique  for  im¬ 
proving  the  quality,  speed,  and  efficiency  of  acquiring  space  surveillance  and  astro¬ 
nomical  imagery.  Finally,  section  2.7  introduces  the  independent  variables  effecting 
target  visibility. 

2.2  Images  and  Objects 

An  image  is  the  result  of  averaging  multiple  short  exposures  (frames)  of  the 
same  object.  The  final  product  of  averaging  together  short  exposure  images  is  re¬ 
ferred  to  as  a  long  exposure  image.  An  object  is  the  actual  source  being  imaged, 
also  referred  to  as  the  target.  Two  types  of  sources  are  point  sources  and  extended 
sources.  A  point  source  is  typically  a  very  bright  single  star  (or  a  laser  glint  off  of  an 
extended  source),  and  is  treated  mathematically  as  a  two-dimensional  delta  function. 
An  extended  source  has  structure  beyond  that  of  a  delta  function.  Examples  of  ex¬ 
tended  sources  are  binary  stars  and  satelHtes.  Figure  2.1  show  the  extended  object 
used  in  this  thesis,  along  with  the  photon  limited  image  before  frame  selection  is 


applied.  Astronomical  image  data  are  essentially  photographs  of  space  objects  such 
as  stars  and  satellites.  This  thesis  deals  primarily  with  point  sources  and  satellites. 

2.3  Atmospheric  Turbulence 

Atmospheric  turbulence  has  been  the  greatest  limiting  factor  for  obtaining 
images  of  space  objects  since  the  invention  of  the  telescope.  Even  Sir  Isaac  Newton 
(1730)  complained  of  its  effects  [17].  With  the  invention  of  photographic  film,  one 
solution  to  the  problem  of  seeing  dim  objects  was  to  expose  the  film  for  a  very  long 
time  and  get  a  brighter  picture  than  could  otherwise  be  seen  with  the  naked  eye. 
But  turbulence  causes  distortion  in  the  images,  as  the  image  moves  and  changes 
shape  during  the  long  exposure.  One  solution  to  this  would  be  to  take  a  very  short 
exposure  and  essentially  freeze  the  turbulence  and  the  motion  of  the  object.  The 
problem  with  short  exposure  images  is  that  they  are  corrupted  by  high  frequency 
modulation.  Consequently,  the  image  has  a  speckled  appearance.  This  problem 
can  be  partially  overcome  by  taking  a  series  of  many  short  exposure  shots  and 
averaging  them  together,  thus  smoothing  out  the  speckles.  These  post  processing 
techniques  are  now  referred  to  as  speckle  imaging.  Therefore,  it  would  seem  that  more 
frames  and  less  exposure  time  is  the  key  to  increasing  the  quality  of  the  final  image. 
However,  exposure  times  are  limited  by  detector  sensitivity  and  and  brightness  levels, 
typically  several  milliseconds.  Therefore,  further  processing  is  required.  Several 
post-processing  techniques  have  been  employed  to  further  increase  the  final  quality 
of  the  image.  The  technique  investigated  here  is  speckle  interferometry  using  frame 
selection,  combined  with  bispectrum  phase  reconstruction. 

2.4  Speckle  Interferometry 

Speckle  interferometry  is  a  method  of  calculating  the  average  power  spectrum 
from  a  collection  image  data  frames.  It  is  a  post-processing  technique  of  combining 
a  series  of  short  exposure  images  and  taking  the  average.  Figure  2.2  illustrates  the 
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Figure  2.1  Extended  Source  Images:  a)  original  image  before  distortion,  b)  im¬ 
age  as  seen  from  a  1  meter  telescope  from  500  kilometers  with  seeing 
conditions  of  To— 7  cm  and  m^=l,  before  frame  selection. 


2-3 


flow  process  of  converting  detected  speckle  image  data  into  the  final  estimate  of  the 
object  [24].  First,  image  data  jfrom  an  unknown  test  object  is  collected  along  with 
reference  data  from  a  nearby  bright  star.  A  laser  glint  off  of  the  object  may  function 
as  the  reference  point  source,  if  a  suitably  near  and  bright  star  is  not  available.  The 
point  source  represents  the  point  spread  function  (PSF)  of  the  test  object,  and  is 
characteristic  only  of  the  atmosphere  and  the  optics.  For  this  reason,  it  is  used  to 
calibrate  for  the  effects  of  the  optics  and  the  atmospheric  distortion.  In  the  absence 
of  the  atmosphere  and  optics,  the  PSFs  of  point  sources  and  an  extended  sources 
axe  pure  delta  functions.  With  an  aberration  free  optical  system,  in  the  absence  of 
atmosphere,  the  PSF  is  an  Airy  disc  [24].  The  size  of  the  disk  is  inversely  proportional 
to  the  size  of  the  entrance  pupil  of  the  optics.  By  introducing  the  atmosphere,  the 
PSF  spreads  even  more.  Aberrations  associated  with  either  the  optical  system,  or  the 
atmosphere,  in  the  form  of  turbulence  cause  the  PSF  to  become  even  more  distorted. 
The  first  step  in  speckle  interferometry  is  to  compute  the  Fourier  transform,  or  the 
image  spectrum,  I{f),  of  each  image,  z(f),  in  the  data  set.  The  spectrum  has  units 
of  inverse  length.  The  Fourier  transform  is  a  mathematical  tool  for  converting  an 
image  to  its  spectrum,  and  is  given  by 

I{f)  =  J i{f)exp{-j2Trf- r}df,  (2.1) 

where  f  is  the  spatial  coordinate,  and  /  is  the  spatial  frequency  coordinate. 

The  inverse  Fourier  transform  is  a  mathematical  tool  for  converting  an  image 
spectrum  back  into  an  image.  Its  form  is 

i(f)  =  J  I{f)  exp{;7r/-  f}  df,  (2.2) 

In  their  two  dimensional  formats,  |r(a;,y)|  =  and  \f{u,v)\  =  +  v^. 

Capital  letters,  such  as  /(/)  and  0[f),  denote  Fourier  transforms  of  a  function. 
These  are  referred  to  here  as  image  and  object  spectra,  respectively.  The  lower  case 
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Figure  2.2  Block  diagram  of  the  speckle  imaging  technique 


equivalents,  z(r)  and  o(f),  refer  to  the  image  and  object,  respectively.  The  result  of 
a  Fourier  transform  is  a  complex- valued  function: 

Hf)  =  exp{i0i(/)}.  (2.3) 

The  term,  |/(/)|,  in  Eq.  (2.3)  is  referred  to  as  the  moduks spectrum  and  exp{j(^j(/)} 
is  the  phase  spectrum  of  the  image  [24]. 

Next,  the  average  modulus  spectrum  is  estimated  by  computing  the  sample- 
based  estimate  of  the  average  of  the  modulus  squared  of  the  spectra  of  the  data 
frames.  The  notations  used  for  the  average  modulus  squared  of  the  Fourier  transform 
of  the  object  and  the  reference  star  are,  respectively,  £^{|/(/)|^}  and  E{\R{f)\'^}, 
where  E{-}  refers  to  the  statistical  expectation  operator  [24]. 

The  image  spectrum  is  related  to  the  object  spectrum,  0{f)  and  the  optical 
transfer  function  (OTF),  'W(/)  by 

lif)  =  0(/)W(/l,  (2.4) 

where  H{f)  is  related  to  the  pupil  function  by  [9] 


.A  f  f-oc-P(^  +  ^,y  +  ^)P(x-^,y-^)dxdy 

ir^P(x,y)dxdy  ’ 

where  di  is  the  distance  between  the  optic  and  the  image  plane,  and  f^  and  fy  are 
the  spatial  frequencies  in  the  x  and  y  directions,  respectively. 

There  is  a  limit  to  how  much  high  frequency  information  from  an  object  can  be 
passed  through  an  optical  system.  This  limit  is  defined  as  the  cutoff  frequency,  /c.  fc 
is  characteristic  of  equipment  and  atmospheric  seeing  conditions.  Spatial  frequency 
corresponds  to  the  sharpness  in  detail  of  an  image  or  object.  For  example,  large  fiat 
surfaces  exhibit  low  spatial  frequencies,  while  ridges  and  sharp  corners  exhibit  high 
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frequencies.  The  lower  the  fc,  the  less  detail  an  image  will  have,  fc  defines  the  point 
in  the  spectrum  where  the  OTF  drops  to  zero.  For  this  reason,  fc  is  also  referred 
to  as  the  OTF  radius.  The  consequence  of  the  frequency  cutoff  is  that  high  spatial 
frequency  information  is  lost.  High  spatial  frequencies  encode  the  fine  details  of  an 
image.  An  image  missing  high  frequency  data  of  an  object  will  appear  whole  but 
blurred.  In  other  words,  lower  frequency  cutoffs  result  in  lower  resolution. 

The  average  modulus  squared  of  the  image  spectrum,  E{\I{f)\‘^},  is  computed 
in  speckle  imaging  using 


=  |0(/)p£{|W(/)p}.  (2.6) 

The  term,  |0(/)p  caabe  pulled  out  of  the  expectation  operator  in  Eq.  {2.6),  since  the 
true  object  spectrum  is  assumed  constant  over  the  duration  of  the  data  collection. 
However,  the  OTF  is  random  due  to  the  random  turbulence-induced  fluctuations 
inherent  in  the  OTF.  Therefore  the  term,  must  remain  inside  the  expecta¬ 

tion  operator.  A  similar  relationship  holds  for  the  reference  star  with  one  further 
simplification,  |0(/)p  =  (7,  a  constant.  For  simphcity,  C  is  set  to  unity  in  Eq.  (2.7). 

£:{|iJ(/)n  =  £{|0(/)W(/)p} 

=  \0{/YE{\H(fY} 

=  E{\H(j')\^}  (2.7) 

The  next  step  is  to  remove  the  effects  due  to  the  average  modulus  squared 
of  the  OTF  from  the  image  power  spectrum.  This  step  is  performed  by  using  the 
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reference  star  data  to  obtain  an  estimate  of  using 

g{|/(/)p}  ^  |o(/)pg{W)P} 

®{|fi(/)p}  Hmfm 

=  |0{/)p.  (2.8) 

Eq.  (2.8)  provides  the  modulus  of  the  object  spectrum.  This  is  the  essence 
of  Labeyrie’s  speckle  interferometry  technique  [24].  In  an  operational  setting,  the 
power  spectrum  of  the  detected  image  is  not  proportional  to  |0(/)p,  but  rather  has 
a  bias  associated  with  noise.  The  actual  modulus  can  only  be  estimated  from  the 
detected  image  by  accounting  for  these  noise  biases.  Noise  bias  effects  on  the  power 
spectrum  axe  discussed  in  the  next  section. 

2.4..I  Noise  Bias  Effects.  Three  sources  of  randomness  are  encountered  in 
obtaining  image  data. 

1.  Turbulence 

2.  Photon  noise 

3.  CCD  read  noise 

Each  of  the  above  sources  has  a  different  statistical  behavior  and  must  be  mathemat¬ 
ically  treated  accordingly.  Turbulence  is  governed  by  turbulence  statistics,  photon 
noise  governed  by  Poisson  statistics,  and  CCD  read  noise  is  governed  by  Gaussian 
statistics.  Turbulence  statistics  are  accounted  for  by  using  Labeyrie’s  technique 
described  in  the  previous  section.  Photon  noise,  also  referred  to  as  shot  noise,  is 
associated  with  the  random  arrival  times  and  locations  of  photo  events  on  the  image 
plan.  When  photon  noise  effects  are  included  in  Eq.  (2.6)  we  find  that  E{|/(/)p}  is 
not  proportional  to  |0(/)p,  but  rather  has  a  bias.  The  model  for  a  detected  image 
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corrupted  with  axiditive  noise  is  given  by  [24] 

(2-9) 

71=1  J)=l 

where  the  first  term  describes  the  photon-limited  image,  and  the  second  term  de¬ 
scribes  the  additive  noise.  The  location  of  the  pixel  in  the  image  is  fp,  and  rip 
is  a  random  variable  representing  the  amount  of  additive  noise  at  the  pixel  [24] . 
The  bias  associated  with  the  additive  noise  is  accounted  for  by  simply  subtracting 
the  total  number  of  photo  events  per  image,  K  from  the  modulus  squared  of  the 
spectrum  of  the  detected  image,  \D{f)\  resulting  in  [3]. 

Qif)  =  \D{f)\'-K,  (2.10) 

where  Q{f)  is  the  unbiased  estimation  of  the  object  spectrum. 

CCD  read  noise  is  a  source  of  noise  which  arises  from  the  readout  electronics 
in  a  CCD  detector.  The  total  number  of  pixels  in  the  image  is  P.  CCD  noise 
is  accounted  for  by  subtracting  Pa^  from  the  unbiased  estimator,  where  is  the 
standard  deviation  of  the  additive  noise  expressed  in  photo  elections  per  pixel  [25]. 
The  equation  for  Q{f)  becomes 

Q{f)  =  \D{f)\''  -K-Pa\  (2.11) 


The  pseudo-Wiener  Filter  [24]  is  an  additional  step  in  the  post-process  of 
speckle  image  data  which  corrects  for  turbulence,  equipment  aberration,  and  filters 
out  some  of  the  high  frequency  noise.  It  is  given  by 


\o{fY  = 


Qobjif) 

Qref{f)  +  Ol/SNRQif)' 


(2.12) 
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The  term  Qrefif)  is  data  from  a  reference  star  and  is  representative  of  the  PSF.  The 
term,  a/SNRQ{f)  reduces  high  frequency  noise.  SNR  is  the  signal-to-noise  ratio, 
described  later  in  section  2.6.2.  The  variable,  a  can  be  adjusted  until  the  desired 
filtering  eflPect  is  achieved. 

The  limitations  to  this  technique  are  that  the  phase  of  the  Fourier  transform 
is  still  required  to  create  an  estimate  of  the  object.  A  phase  estimator  method  such 
as  the  bispectrum  technique,  must  be  employed  to  estimate  the  phase  of  the  Fourier 
transform  of  o(r). 

2.5  Bispectrum  Calculations  . 

Since  speckle  interferometry  provides  only  the  modulus  of  the  object  spectrum, 
|0(/)|,  bispectrum  phase  reconstruction,  reconstructs  the  phase  spectrum  of  the 
object,  4>o{f),  [24].  The  bispectrum  is  a  specialized  moment  of  the  image  spectrum 
which  contains  high  spatial  frequency  phase  information. 

2.5.1  Bispectrum  Definition.  The  bispectrum  of  an  image  is  defined  as  [13] 

A)  =  +  h),  (2.13) 

where  /i  and  /2  are  any  two  points  in  the  image  spectrum.  The  phase  of  the  bis¬ 
pectrum  encodes  the  phase  of  the  object  as  point-to-point  phase  differences  [24]. 
The  object  phase  spectrum,  is  obtained  from  sample  based  estimates  of  the 

bispectrum.  Note,  that  since  /i  and  /2  are  functions  of  spatial  x  and  y  coordinates, 
the  bispectrum  is  a  four  dimensional  data  object.  Hermite  symmetry  of  the  image 
spectrum  results  in  an  eight-fold  bispectrum  symmetry  expressed  by  [13] 


B(/l,  /j)  =  Bih,  /i)  =  B(/,  -  4  /I)  =  (2.14) 
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The  object  phase  spectrum  is  encoded  in  the  phases  of  the  bispectrum  according  to 
the  following  relationship  [24]: 

B(h,h)  =  \o(h)\\o(!2)\\o(h-^h)\ 
m(h)\\H(h)\\H(h  +  h)\ 

X  exp{jl4>4fi)  +  +  A)  +  <Pn(fi) 

+<pH{f2)  —  <hi{h  + h)]}>  (2.15) 

where  the  exp{j[^<,(A)  +  4>o{f2)  —  <Po{fi  +  ^)]}  terms  contain  the  encoded  object 
phase  spectrum  [24].  Any  corruption  in  the  phase  spectrum  by  the  random  phase 
spectrum  of  the  OTF  is  overcome  by  taking  the  average  of  several  data  points.  The 
object  spectrum  is  then  related  to  the  expected  value  of  the  bispectrum  by  [24] 

/!)}  =  0(/i)0(/2)0'(/i  +  /2)B{«(/i)W(/j)«-{/i  +  h)},  (2.16) 

where  the  phase  of  E{H{fi)H{f2)'H*{fi  +  ^)}  is  zero. 

The  symmetric  nature  of  the  bispectrum  helps  reduce  the  required  number 
of  calculations.  Roggemann  points  out  the  importance  of  reducing  the  size  of  the 
bispectrum  data  object  [24].  After  considering  all  the  possible  combinations  of  /i 
and  /2  and  the  eight-fold  symmetry  of  the  bispectrum,  a  128  x  128  pixel  image  results 
in  128^  X  (128^  —  l)/8  =  33,032,192  calculations.  One  solution  is  to  constrain  the 
range  of  values  {fi  +  are  allowed  to  take  [24].  This  also  argues  a  case  for  the 
benefits  of  frame  selection  as  a  means  to  data  reduction. 

S.5.2  Insensitivity  to  Tilt  Error.  The  bispectrum  technique  has  an  advan¬ 
tage  over  an  alternate  phase  estimation  technique,  known  as  the  Knox-Thompson 
cross  spectrum  technique  [26],  in  that  the  bispectrum  is  insensitive  to  random  mo¬ 
tion  of  the  image  centroid  [24].  The  necessity  to  correct  for  tilt  by  recentering  all  the 
images  is  eliminated.  Any  tilt  errors,  or  linear  shifts  in  the  centroid  appear  as  linear 
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phase  shifts  in  the  spectrum.  For  example,  if  a  shift  in  the  image  is  represented 
by  i(r  +  Ar^,  the  spectrum  becomes  exp{— j27r/*-  Af}I{f).  The  bispectrum  of  the 
shifted  image  is  given  by  [24] 

Bs{flj2)  =  I{fl)l{f2)l%fl  +  f2) 

X  exp{-j27r[/i  •  Af  +  ^  •  Af-  (/i  +  ^)  •  Ar]} 

=  (2.17) 

Therefore,  the  bispectrum  of  the  shifted  image  is  equivalent  to  the  bispectrum  of  the 
unshifted  image. 

2.5.3  Noise  Bias  Effects.  The  unbiased  estimator  for  the  bispectrum  is 
obtained  from  the  spectrum  of  the  detected  image,  D{f),  and  is  corrected  for  additive 
noise  effects.  Its  equation  is  given  by  [24] 

=  B(/i)D(A)D*(/i  + /2)  -  |Z)(/i)P 

-|D(A)P  -  \D(fi  +  hf  +  2K  +  3F<T^  (2.18) 


where  K,  P,  and  a  are  as  defined  in  section  2.4.1. 


2.5.4  Phase  Error.  The  estimated  phase  spectrum,  4>o{f),  is  derived  from 
the  bispectrum.  The  specific  method  of  reconstructing  the  phase  spectrum  of  the 
object  from  the  bispectrum  is  detailed  in  section  3.4.1.  As  a  measure  of  performance, 
(j>o{f)  is  compared  to  the  actual  phase  spectrum,  </>o{f),  of  the  true  object  in  Fig.  2.1 
by 


Mff  =  (Uf)  -  Mf)?,  (2.19) 


where  0«(/)  is  the  phase  error  in  units  of  radians  squared.  The  phase  spectrum  is 
then  combined  with  the  power  spectrum  and  inverse  Fourier  transformed  into  the 
image. 
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The  relative  phase  error  performance  of  a  case  with  frame  selection  is  compared 
to  that  without  by 

^a(/)  =  4‘eif)  100%  ~  (2.20) 

where  (pe{f)ioo%  represents  the  phase  error  of  the  100%  case,  and  <f)e{f)FS  represents  a 
frame  selected  case.  Areas  in  the  spectrum  where  performance  is  improved  as  a  result 
of  frame  selection,  will  have  positive  4>a  values.  Conversely,  where  the  performance 
declines  as  a  result  of  frame  selection,  4>a  values  are  negative. 

2.6  Frame  Selection 

This  section  discusses  the  frame  selection  process.  It  is  often  observed  that 
for  any  given  set  of  image  data,  there  are  some  images  that  are  simply  better  than 
others  [4,  5,  7,  18,  22,  23,  25].  The  goal  of  frame  selection  is  to  devise  a  method 
for  sifting  through  large  sets  of  image  data  frames  and  selecting  the  fewest  number 
of  the  best  frames  that  maintains  or  improves  the  quality  of  the  final  image.  Frame 
selection  is  a  three-step  process  for  choosing  a  subset  of  images  for  further  processing. 

1.  Calculate  a  frame  quality  metric  value  for  each  frame  in  the  data  set. 

2.  Sort  the  images  in  the  data  set  from  highest  to  lowest  based  on  the  quality 
metric  value. 

3.  Select  the  highest  value  images  for  processing. 

Captains  Craig  Stoudt  and  Steve  Ford  have  done  extensive  research  in  the  areas 
of  frame  selection  combined  with  adaptive  optics  [5,  25].  They  found  that  under 
favorable  brightness  levels,  frame  selection  improves  adaptive  optic  performance. 

An  effective  implementation  of  frame  selection  depends  on  the  size  of  the  data 
base  to  be  processed.  For  extremely  large  data  bases,  a  computer  frame  selection 
algorithm  containing  a  simple  but  effective  quality  metric,  may  be  used  to  effectively 
reduce  the  amount  of  data.  The  selected  subset  may  then  be  processed  using  a 
bispectrum  phase  reconstruction  algorithm. 
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The  basis  of  frame  selection  is  to  choose  the  best  quality  or  highest  resolution 
frames  for  averaging  such  that 

1  M 

W)  =  ]gEWl,  (2-21) 

where  Is{f)  is  the  average  image  frequency  spectrum  due  to  frame  selection,  Ii{f)  is 
the  image  frequency  spectrum  of  an  individual  frame  in  the  frame  selection  subset, 
M  is  the  number  of  frames  in  the  selected  subset,  and  /  is  a  spatial  frequency  [5]. 
Figure  2.3  provides  a  block  diagram  of  the  process  associated  with  Eq.  (2.21)  [23]. 
N  images  are  collected  and  Fourier  transformed.  The  images  are  then  sequentially 
rank-ordered  from  0  to  1,  according  to  a  quality  metric  determination  and  the  total 
number  of  frames.  The  best  image  is  assigned  the  number  1,  while  the  worst  image 
receives  0.  The  images  that  are  above  a  user  defined  cutoff  (between  0  and  1)  are 
maintained  for  further  processing  and  those  below  are  discarded.  The  user  defines 
the  cutoff  through  the  frame  selection  rate  (FSR).  This  is  discussed  in  the  next 
section. 

2.6.1  Frame  Selection  Rate  (FSR).  The  frame  selection  rate,  FSR,  is  the 
ratio  of  the  number  of  frames  M  in  the  selected  subset  to  the  number  of  frames  N 
in  the  original  ensemble.  Therefore, 


FSR  =  — ,  M<N.  (2.22) 

The  FSR  is  a  percentage  value.  An  FSR  of  70%  means  the  selected  subset  consists 
of  the  best  70%  of  the  frames  from  the  original  data  set  based  on  the  frame  quality 
metric  values  [5].  In  the  rank  order  system,  ail  frames  ranked  below  (1-FSR) — in 
the  70%  case,  0.299  and  below — are  discarded. 

2.6.2  Signal-to-Noise  Ratio  (SNR).  As  an  aid  to  the  user,  a  metric  for 
judging  the  quality  of  the  power  spectrum  estimate,  the  signal-to-noise  ratio  (SNR), 
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Figure  2.3  Block  diagram  of  the  frame  selection  technique 


is  calculated  to  compare  the  image  quality  of  the  total  set  versus  the  selected  sub¬ 
set  [5].  The  SNR  of  the  image  spectrum  is  defined  by 


SNRoif)  = 


EjQif)} 


(2.23) 


The  SNR  of  a  frame  selected  case  is  compared  to  that  of  the  100%  (without 
frame  selection)  case  by 


SNRgain  = 


SNRfs 
S  NRiqo%  ’ 


(2.24) 


where  SNRgain  is  the  SNR  gained  as  a  result  of  frame  selection,  SNRps  is  the  SNR 
resulting  from  frame  selection,  and  SNRiqo%  is  the  SNR  resulting  from  not  using 
frame  selection.  Areas  in  the  spectrum  where  performance  is  improved  as  a  result 
of  frame  selection,  will  have  SNRgain  values  greater  than  1.  Conversely,  where  the 
performance  declines  as  a  result  of  frame  selection,  the  SNRgain  values  less  than  1. 


One  unavoidable  consequence  of  frame  selection  is  that  it  changes  the  statis¬ 
tics  of  the  data  set.  Frame  selection  has  been  shown  to  increase  the  total  SNR  in 
some,  but  not  all  circumstances  [5,  25].  This  is  because,  although  frame  selection 
is  designed  to  improve  the  data  set  by  discarding  the  worst  frames,  there  comes  a 
point  when  throwing  out  more  bad  data  means  throwing  out  valuable  information 
as  well.  This  introduces  an  important  limitation  explained  by  the  statistical  nature 
of  noise  reduction  through  averaging  [2].  The  SNR  of  a  set,  SNRtotah  containing  N 
data  points  has  been  shown  to  increase  by  a  factor  y/N  over  that  of  a  single  data 
point  [2].  The  average  SNR  of  a  data  set,  S N Ravg{totai),  is  representative  of  a  single 
data  point.  This  relationship  is  given  by 


SNRtatal  =  ^SNRangitotal)  (2.25) 


2-16 


In  the  case  of  frame  selection,  where  M  frames  are  retained,  M  <  N,  the  SNR 
becomes 

SNR^y,f,fet  =  ''/MSNR  avg(subset)  (2.26) 

Therefore,  in  order  for  image  quality  to  benefit  from  frame  selection,  the  SNR  of 
the  subset  must  improve  by  at  least  a  fax;tor  of  However,  if  the  primary  goal 
in  frame  selection  is  to  reduce  data,  the  total  SNR  may  be  allowed  to  decrease 
slightly.  After  all,  the  ultimate  test  of  image  quality  is  by  visual  inspection.  Slight 
changes  in  SNR  axe  often  unnoticeable  to  the  human  eye.  In  this  situation,  the 
restriction  may  be  relaxed  to  a  user  defined  tolerance  level.  is  the  minimum 
SNR  reduction,  for  frame  selection  to  be  more  effective  than  selecting  the  frames 
at  random,  or  just  simply  collecting  a  smaller  set  of  image  data  in  the  first  place. 
Fig.  2.4  is  a  typical  SNR  plot  of  a  point  source  showing  the  resulting  SNR  plots  of 
data  that  has  been  randomly  selected  without  the  benefit  of  a  sharpness  metric.  It 
shows  that  the  gradual  reduction  in  the  SNR  as  the  number  of  randomly  selected 
data  is  reduced. 

The  limitation  to  the  SNR  metric,  is  that  it  measures  the  composite  image 
based  only  on  the  power  spectrum.  An  additional  metric  is  required  to  measure 
the  composite  image  based  on  the  phase  spectrum.  This  is  the  phase  error  metric, 
discussed  in  section  2.5.4. 

There  are  several  methods  available  for  sifting  through  large  amounts  of  data 
for  frame  selection.  One  can  computerize  the  process,  after  obtaining  an  effective 
quality  metric.  The  next  section  discusses  some  specific  quality  metrics. 

2.6.3  Quality  Sharpness  Metrics  for  Individual  Frames.  A  quality  metric 
is  an  algorithm  which  results  in  a  single  numerical  value  when  applied  to  an  image. 
The  magnitude  of  a  quality  metric  is  directly  related  to  a  designated  quality  or 
feature  of  the  image,  i{f),  obtained  in  Eq.  (2.2).  These  metrics  are  known  to  take 
their  maximum  value  for  a  diffraction  limited  image.  The  quality  metrics  used  in 
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Figure  2.4  Example  SNR  plot  illustrating  the  theoretical  relationship  of  SNR  to 
randomly  selected  data,  mj,  =  1,  ro  =  12  cm. 
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this  thesis  are  derived  from  Muller  and  Buffington’s  eight  sharpness  functions.  The 
complete  list  and  a  brief  discussion  is  given  below: 


Si  =  J  j  if)  dr 

(2.27) 

5'2  =  i(0,0) 

(2.28) 

S'a  =  y  y  m{f)i{f)  dr 

(2.29) 

S,=  f  f  '  df 

J  J  dx'^dy^ 

(2.30) 

Sz  -  J  j  i"(r)df 

(2.31) 

Sq  =  -  J  J  i(f)r^  dr, 

(2.32) 

87  =  —  J  J  fn[i(r)]z(r)  dr 

(2.33) 

'S  =  -  J  j  K(f^)  -  o(f)p  dr. 

(2.34) 

is  the  maximum  structural  content  metric.  m(f)  is  known  as  the  mask  function. 
If  m(f)  is  a  good  replica  of  the  undistorted  image,  S^,  Eq.  (2.29),  is  a  good  sharpness 
definition  which  reduces  to  ^i,  Eq.  (2.27)  when  distortion  is  removed  [16].  S2  works 
well  only  for  bright  objects.  S4  is  the  total  differential  metric.  87  is  known  as  the 
maximum  entropy  metric.  Ss  is  the  fidelity  defect  function,  which  requires  prior 
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knowledge  of  the  object,  o(r).  Sg  is  only  useful  in  an  experimental  setting  as  a 
comparison  check  for  judging  the  performance  of  other  metrics.  These  metrics  reach 
maximum  values  for  unaberrated  images  [16].  In  the  form  of  computer  algorithms, 
along  with  the  associated  maximum  values,  these  metrics  are  tools  for  comparing 
the  results  of  test  data  and  assigning  numerical  values  for  determining  the  relative 
quality  of  an  image  within  a  given  set  of  image  frames.  The  actual  metrics  used  in 
this  thesis  are  discussed  in  section  3.3.2.  The  next  section  discusses  parameters  that 
may  vary  the  overall  quahty  of  speckle  images. 

2. 7  Independent  Parameters 

In  this  section,  two  key  independent  parameters  in  speckle  imaging,  target 
brightness  levels  and  turbulence  are  discussed.  These  are  referred  to  as  visual  mag¬ 
nitude  and  Fried’s  parameter  of  seeing  conditions,  respectively. 

2.7.1  Visual  Magnitude.  Visual  magnitude,  is  a  system  employed 
by  astronomers  to  describe  the  relative  brightness  of  objects  in  a  night  sky  [25]. 
The  average  number  of  photons,  is  calculated  based  on  this  number,  which  is 
dependent  on  the  rate  of  incident  photons  per  unit  area,  or  photon  flux,  P;  the  total 
area,  or  the  telescope  pupil  size,  Apup',  and  the  integration  time  of  the  processor,  r. 

Terry  Gray  of  the  Phillips  Laboratory,  Kirtland  AFB,  NM,  wrote  a  program 
called  VMAG  for  producing  values  for  P  based  on  a  given  target  wavelength,  band¬ 
width,  and  telescope  transmission  factor  [10].  These  values  are  for  a  standard  atmo¬ 
spheric  transmission  factor,  «  of  1.  Therefore,  the  simulation  must  account  for  less 
than  perfect  atmospheric  conditions,  by  including  k  in  the  equation  as  follows: 

K  =  {Ap,p){P){t){k)  (2.35) 

Table  2.1  lists  visual  magnitudes  and  corresponding  photon  intensities  of  typi¬ 
cal  night  time  objects  for  a  1  meter  aperture,  r  of  1.8  ms,  and  a  wavelength  of  500  nm 
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Table  2.1  Visual  magnitude  and  photo  events  for  common  sky  objects 


OBJECT 

ruj, 

K 

Sun 

-26.7 

1.3 

lO^o 

Full  Moon 

-12.5 

4.7 

10^ 

Venus 

-4.4 

2.35 

10^ 

Jupiter 

-2.5 

4.11 

10® 

Sirius 

-1.5 

1.59 

10® 

Artificial  Satellite  {typical) 

-0.96 

1.00 

10® 

Polaris 

+2.2 

5.54 

10^ 

with  a  10%  bandwidth  [5].  Note  that  the  objects  get  brighter  as  decreases  in 
value. 


2.7.2  Pried  Seeing  Parameter.  As  ground-based  astronomical  telescopes 
become  increasingly  larger  and  the  optics  become  more  powerful  and  less  flawed, 
the  ability  to  resolve  distant  images,  such  as  satellites,  stars,  and  neighboring  galax¬ 
ies,  becomes  limited  not  by  diffraction  nor  aberration  due  to  the  optics,  but  by 
atmospheric  distortion,  due  to  turbulence  [16].  The  Fried’s  seeing  parameter  ro,  is 
a  mathematical  measure  of  this  limit.  For  a  given  atmospheric  condition,  To,  the 
atmospheric  coherence  diameter,  is  a  figure  of  merit  for  the  ability  to  resolve  images 
through  it.  For  large  telescopes  this  limits  the  resolution  to  that  obtainable  by  a 
telescope  of  diameter  r^  [22].  The  probability  of  obtaining  diffraction  hmited  images 
is  exponentially  dependent  on  the  ratio  squared  of  the  pupil  diameter,  D,  divided 
by  ro,  as  shown  by 

(2.36) 

As  the  diameter  of  the  optic  becomes  much  greater  than  the  seeing  parameter  then 
the  image  is  said  to  be  limited  by  atmospheric  distortion.  If  D  •<  ^o,  fc  =  DfXdi, 
If  D  ^  ro,  fc  =  folXdi,  where  D  is  the  telescope  diameter,  fc  is  the  OTF  radius, 
discussed  in  section  2.4,  A  is  the  mean  wavelength  of  the  detected  image,  and  dj  is 
the  distance  between  the  aperture  and  the  image  plane.  Typical  values  of  ro  at  good 
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observatories  is  approximately  10-30  cm  [24],  at  visible  wavelengths.  The  current 
state-of-the-art  technology  far  surpasses  this  in  the  ability  to  create  optical  systems. 
But  the  Air  Force’s  desire  to  have  close  up  views  of  every  satellite  in  orbit  would 
require  far  more  resources  than  is  currently  feasible.  Therefore,  an  alternate  solution 
to  this  problem  is  to  compensate  for  turbulent  distortion  through  adaptive  optics 
or  to  deconvolve  the  images  through  the  use  of  Fourier  optics,  in  particular,  speckle 
interferometry,  combined  with  a  systematic  process  of  frame-selection.  The  next 
section  discusses  some  recent  AFIT  accomplishments  in  the  application  of  frame 
selection  to  adaptive  optics. 

2.8  Recent  Developments  in  Frame  Selection 

In  1993,  Craig  Stoudt,  an  AFIT  graduate  student,  investigated  the  perfor¬ 
mance  of  frame  selection  by  simulating  a  1  meter  fully  compensated  adaptive  optics 
telescope  [25].  He  evaluated  three  of  Muller  and  Buffington’s  original  quality  metrics, 
Si,  Si,  and  Ss,  as  well  as  a  new  metric  Ssi  defined  as 

Ssi  =  I  JSoTriDWYdf  (2.37) 

where  SoTpif)  is  the  support  of  the  OTF.  Since  anything  out  side  of  the  OTF  radius 
is  pure  noise,  this  term  defines  limits  to  the  integration  so  that  everything  outside 
the  radius  is  ignored.  Stoudt’s  use  of  computer  simulation  allowed  a  much  broader 
examination  of  the  effects  of  key  independent  parameters  than  those  conducted  by 
earlier  researchers.  Image  spectrum  SNR  was  used  as  the  primary  performance 
measure.  He  showed  that 

1.  Frame  selection  rates  of  60-70%  improves  the  image  spectrum  SNR  by  ten  to 
fifteen  percent  at  mid  spatial  frequencies  for  both  point  and  extended  sources. 

2.  Lower  selection  rates  (10-20%)  result  in  sharper  but  noisier  images. 
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3.  For  equivalent  seeing  conditions  and  light  levels,  frame  selection  is  of  greater 
benefit  to  point  sources  than  extended  objects. 

4.  Frame  selection  effectively  corrects  for  poor  seeing  conditions,  but  is  not  as 
effective  in  correcting  for  photon  noise  and  CCD  camera  noise. 

5.  The  Ssi  quality  metric  produced  the  highest  overall  image  spectrum  SNR  gain 
compared  to  the  other  metrics. 

In  1994,  Steve  Ford,  also  an  AFIT  graduate  student,  demonstrated  frame  se¬ 
lection  on  actual  Air  Force  Maui  Optical  Site  (AMOS)  imagery.  He  has  shown  that: 

1.  Frame  selection  provides  improved  image  spectrum  SNR  by  discarding  the 
worst  manifestations  of  atmospheric  turbulence.  This  results  in  reduced  image 
spectrum  variance  and  higher  SNR. 

2.  For  point  source  objects,  the  minimum  visual  magnitude  for  SNR  gain  due 
to  frame  selection  was  m^,  =  +7  —  8  depending  on  the  adaptive  optics  (AO) 
system  modeled.  An  example  extended  object,  a  satellite  model,  required  a 
much  brighter  minimum  visual  magnitude  of  —  +A.  In  general,  average 
seeing  conditions,  represented  by  tq  =  10  centimeters,  benefited  the  most  from 
frame  selection  processing.  This  increased  SNR  gain  is  due,  in  part,  to  larger 
OTF  variance  for  average  seeing  conditions  when  compared  to  poor  or  very 
good  conditions. 

3.  CCD  camera  noise  effects  play  a  significant  role  in  estabhshing  frame  selection 
performance  limits. 

4.  When  objects  are  sufficiently  dim,  photon  and  CCD  camera  noise  effects  dom¬ 
inate  the  single  frame  SNR  expression.  All  available  frames  in  the  data  set 
are  needed  to  reduce  this  effect  through  averaging.  Therefore,  frame  selection 
cannot  provide  image  improvement. 

5.  Point  spread  functions  (PSFs)  are  nearly  identical  for  point  source  and  ex¬ 
tended  objects  after  frame  selection  processing.  This  fact  allows  the  optimal 
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use  of  deconvolution  techniques,  such  as  inverse  and  Weiner  filters,  to  sharpen 
images  after  frame  selection. 

6.  Short  exposure  times  are  necessary  to  get  the  full,  predictable  benefits  of  frame 
selection  processing  on  real  imagery.  Improvements  in  SNR  gain  were  demon¬ 
strated  for  reference  star  and  satellite  objects  collected  by  the  AMOS  facility. 
Experiments  using  this  data  were  not  readily  repeatable  due  to  the  60-99  mil¬ 
lisecond  exposure  times  used  for  image  collection. 

2.9  Limits  of  Current  Knowledge 

The  Air  Force  uses  adaptive  optics  (AO),  to  compensate  for  the  random  phase 
aberrations  due  to  atmospheric  turbulence.  Unfortunately,  AO  cannot  completely 
compensate  for  aU  phase  fluctuations  due  to  the  atmosphere  [27].  Therefore,  some 
form  of  post-detection  processing  is  needed.  Recent  research  at  the  Air  Force  Insti¬ 
tute  of  Technology  (AFIT)  has  shown  that  images  from  AO  systems  can  be  improved 
using  frame  selection  [25,  5].  However,  in  many  situations  adaptive  optics  are  cost 
prohibitive.  Therefore,  an  effective  speckle  imaging  technique  using  frame  selection 
is  desirable.  In  the  past,  researchers  have  improved  images  from  passive  optics  sys¬ 
tems  using  frame  selection  [4,  7,  18].  However,  so  far,  nothing  is  known  about  frame 
selection  as  it  applies  to  speckle  imaging. 

2.10  Summary 

This  chapter  laid  the  foundation  for  understanding  the  problem  of  processing 
Air  Force  space  surveillance  and  astronomical  image  data.  Atmospheric  turbulence 
imposes  a  fundamental  limit  on  optical  system  performance.  Speckle  imaging  tech¬ 
niques  enhance  image  quality  to  some  degree,  but  requires  the  post-detection  process¬ 
ing  of  several  hundred  image  data  frames  to  create  a  singe  object.  Frame  selection 
is  a  promising  post-detection  processing  technique  that  offers  dramatic  reduction 
in  data  without  significant  loss  in  image  spectrum  SNR,  or  increase  in  bispectrum 
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phase  error.  A  complete  investigation  of  frame  selection  performance  over  a  substan¬ 
tial  range  of  performance  conditions  and  object  brightness  levels  is  presented  in  this 
thesis.  This  effort  includes  demonstrating  the  quality  of  images  after  frame  selection 
and  sharpening  with  deconvolution  techniques.  This  chapter  provided  the  basic  ter¬ 
minology  used  in  frame  selection,  speckle  imaging  and  the  bispectrum  technique.  It 
outhned  the  performance  metrics,  SNR  and  bispectrum  phase  error,  used  to  evalu¬ 
ate  the  quality  of  composite  images,  and  defined  the  independent  parameters  that 
are  varied  in  the  experiments.  The  methodology  for  conducting  this  investigation  is 
developed  in  the  next  chapter. 
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III.  Approach  and  Methodology 


3.1  Introduction 

Chapter  two  provided  a  brief  overview  of  the  atmospheric  imaging  problem  and 
several  possible  solutions.  A  methodology  for  studying  the  performance  of  frame 
selection  and  bispectrum  phase  reconstruction,  as  they  pertain  to  uncompensated 
speckle  images,  is  developed  in  this  chapter.  Section  3.2  outlines  the  steps  in  the 
computer  simulation  for  both  frame  selection  and  bispectrum  phase  reconstruction. 
Section  3.3  identifies  the  independent  variables  used  in  these  experiments,  including 
the  various  metrics  employed  and  the  justification  for  the  use  of  the  Ssi  frame  quality 
metric. 

3.2  Computer  Simulation 

Computer  simulation  is  a  critical  tool  for  simulating  speckle  imaging,  study¬ 
ing  the  performance  of  frame  selection,  and  investigating  the  bispectrum  technique. 
The  user  input  data  file  allows  convenient  manipulation  of  key  independent  param¬ 
eters.  The  simulation  provides  useful  image  spectrum  statistics,  such  as  the  image 
spectrum  SNR,  and  bispectrum  phase  error,  which  are  measures  of  performance. 
The  imaging  and  frame  selection  code  used  in  this  thesis  was  developed  by  Michael 
Roggemann  [20]  and  further  modified  by  Craig  Stoudt  [25].  The  bispectrum  code 
was  also  developed  by  Michael  Roggemann  [20]. 

Fig  3.1  is  a  block  diagram  of  the  computer  simulation  of  the  speckle  imaging 
process.  First,  the  parameters  are  read  into  the  program  by  a  user  generated  data  file. 
Next,  the  object,  o(r),  is  created,  scaled,  and  Fourier  transformed  into  the  spectrum, 
0{f).  For  purposes  of  this  investigation — simulating  speckle  imaging,  under  ideal 
equipment  conditions — adaptive  optics  was  turned  off  and  the  CCD  read  noise  was 
set  to  zero.  A  series  of  short  exposure  images  are  created  by  a  five  step  process: 


3-1 


•  A  random  phase  screen  modeling  atmospheric  turbulence  is  created. 

•  The  pupil  is  imposed  to  create  the  diffraction  effects. 

•  The  OTF,  or  H{f)  is  calculated. 

•  The  image  spectrum,  /(/),  is  computed  by  multiplying  the  OTF  with  the 
object  spectrum,  and  a  photon-limited  image,  i(f),  corresponding  to  I{f)  is 
created. 

•  The  image,  z{f),  is  Fourier  transformed  to  obtain  the  detected  image  spectrum, 
£>(/),  and  the  total  number  of  photo  events  in  the  image,  K,  is  calculated. 

Once  N  images  are  created,  the  frame  selection  algorithm  sorts  the  images  based  on 
quality  and  selects  M  images  to  be  retained  for  further  processing.  Each  of  the  M 
images  are  sent  through  a  two  step  process: 

•  The  unbiased  power  spectrum,  Q{f),  and  bispectrum,  Bu{fi,f2)^  are  com¬ 
puted. 

•  The  statistics  of  Q{f)  and  are  accumulated. 

Finally,  the  average  unbiased  power  spectrum  and  bispectrum  are  calculated  based 
on  the  subset.  Prom  here  the  phase  spectrum  is  calculated,  compared  to  the  phase 
spectrum  of  the  original  object,  0(/),  for  error  analysis,  and  then  the  image  is 
reconstructed. 

3.3  General  Parameters 

For  comparison,  image  data  statistics  and  the  bispectrum  are  accumulated  in 
two  stages:  once  before  and  once  after  frame  selection.  This  data  is  used  to  determine 
SNR,  SNR  gain,  phase  error,  and  comparative  phase  error  difference  results  of  the 
composite  image  both  with  and  without  the  benefits  of  frame  selection.  SNR  is 
discussed  in  section  2.6.2,  and  phase  error  is  discussed  in  section  2.5.4.  Initially,  the 
program  creates  the  image  composed  of  a  series  of  data  frames,  then  accumulates 
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Figure  3.1  Block  diagram  of  the  speckle  imaging  computer  simulation 
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statistics,  calculates  the  bispectrum,  and  computes  the  quality  metric  values  for  each 
short  exposure.  These  quality  metric  values  are  stored  in  a  file.  Next,  the  program 
reads  the  quality  metric  file  and  sorts  these  values  from  highest  to  lowest.  The  frame 
selection  routines  continue  by  calculating  a  minimum  quality  metric  value  based  on 
the  FSR  indicated  by  the  user.  Then,  the  program  recreates  only  the  short  exposure 
images  with  quality  metric  value  greater  than  or  equal  to  the  minimum  value.  Next, 
the  program  accumulates  a  new  set  of  statistics  and  a  new  mean  bispectrum  based 
on  the  selected  subset.  Finally,  the  results  of  the  simulation  are  deconvolved  and 
measured  against  the  original,  unaberrated,  scaled  object  in  terms  of  phase  error. 
The  phase  errors  are  radially  averaged  as  a  function  of  spatial  frequency,  normalized 
to  the  OTF  radius  defined  by  the  fc  as  was  discussed  in  section  2.4. 

Four  different  parameters  were  examined  to  determine  the  behaAdor  of  frame 
selection  for  speckle  imaging.  Two  of  these  are  user  defined:  the  effectiveness  of 
various  frame  selection  rates;  and  five  different  quality  metrics.  The  other  two  are 
independent  variables  describing  brightness  levels,  and  seeing  conditions.  Other 
parameters  remained  fixed  throughout  the  investigation.  These  include  the  following: 

•  number  of  short  exposure  images  =  300 

•  optic  diameter  =  1  m 

•  CCD  noise  level  =  0 

•  target  range  =  500  km 

•  target  size  =  10  square  meters 

•  mean  wavelength  =  700  nm 

•  bandwidth  =  10% 

•  integration  time  =  10  ms 

•  atmospheric  transmittance  =  10% 

•  telescope  transmission  factor  =  10% 
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3.3.1  Number  of  Short  Exposure  Images.  The  total  number  of  frames  used 
for  each  simulation  run  was  300.  This  number  is  sufficiently  high  enough  to  provide 
an  adequate  sample  base  from  which  statistical  data  can  be  drawn.  The  importance 
of  having  a  large  data  base  is  emphasized  in  section  2.6.2.  Yet,  it  is  small  enough 
to  permit  the  rapid  completion  of  a  series  of  computer  runs  within  reasonable  time. 
Another  point  to  consider  is  that  the  number  of  images  in  a  data  set  is  limited  by 
the  motion  of  the  object  as  it  moves  across  the  sky.  A  satellite  may  pass  through 
the  field-of-regard  of  a  ground-based  telescope  in  a  matter  of  minutes,  limiting  the 
number  of  frames  collected  to  a  few  hundred  [5]. 

3.3.2  Quality  Sharpness  Metrics.  Frames  were  selected  using  the  image 
sharpness  metrics  developed  by  Muller  and  Buffington  [16].  Refer  to  section  2.6.3  for 
a  complete  list  and  a  mathematical  description.  Recall,  from  Eq.  (2.37)  is  the 
modified  version  of  one  of  the  original  eight,  created  at  AFIT  by  Capt  Craig  Stoudt 
in  1993  [25].  This  thesis  examines  the  following  sharpness  metrics: 

•  Structural  Content  (5i  ) 

•  Partial  Differential  (54) 

•  Maximum  Entropy  (^7) 

•  Maximum  Spectrum(5’5i)  (a  modified  ^i) 

•  Fidelity  Defect  {S^) 

The  5^1  metric  is  introduced  in  section  2.8.  Its  equation  is  repeated  here: 

Ssi  =  j  JsoMf)wYdf  (3.1) 

This  is  equivalent  to  the  Si  metric,  except  the  Sqtf  term  limits  the  integration  to 
within  the  OTF  radius.  Stoudt  concluded  that  Ssi  was  superior  in  the  widest  set 
of  circumstances  when  combined  with  adaptive  optics  [25].  Ford’s  investigation  was 


3-5 


narrowed  to  the  performance  limitations  of  this  specific  metric  [5].  However,  this 
investigation  examines  the  other  metrics  to  compare  their  performance  on  speckle 
images. 

3.3.3  Frame  Selection  Rate  (FSR).  Five  different  FSRs  were  tested:  100%, 
80%,  70%,  60%,  and  50%  FSR.  The  FSR  can  be  selected  by  entering  the  desired 
rate  in  the  input  data  file.  The  100%  case  represents  the  baseline  against  which  all 
other  FSRs  are  compared.  Previous  works  have  shown  that  frame  selection  combined 
with  adaptive  optics  has  optimal  values  dependent  upon  the  selected  values  of  the 
independent  variables.  Optimal  FSRs  ranged  from  60-80%  [5,  25].  Therefore,  it  was 
the  initial  intent  of  this  thesis  effort  to  narrow  the  investigation  down  to  the  60%,  70% 
and  80%  cases.  However,  results  from  the  speckle  imaging  data  has  indicated  very 
little  change  in  performance  in  any  of  the  FSRs  tested.  Therefore,  the  investigation 
was  extended  to  include  the  50%  FSR  as  well. 

3.3.4  Visual  Magnitude.  This  investigation  explores  objects  with  values 
of  +1,  +3,  +5,  and  +7,  which  are  all  dimmer  than  most  items  listed  in  table  2.1. 
The  visual  magnitude  of  an  object  can  be  selected  in  the  simulation  by  changing  the 
value  of  D  and  providing  a  value  of  P  obtained  from  a  table  of  flux  values  [10].  Both 
Ford  and  Stoudt  conclude  that  out  of  all  the  independent  parameters,  mi,  has  the 
greatest  effect  on  frame  selection  performance  [5,  25]. 

3.3.5  Fried  Seeing  Parameter.  The  primary  independent  variable  related  to 
the  mean  and  the  variance  of  the  OTF,  is  the  relative  seeing  conditions,  characterized 
by  the  parameter  ro  [25].  The  three  values  tested  in  this  thesis  represent  a  wide 
range  of  operating  conditions: 

•  7  cm  (poor  seeing  conditions) 

•  12  cm  (average  seeing  conditions) 

•  17  cm  (excellent  seeing  conditions). 
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The  poor  seeing  conditions  cause  the  greatest  variance  in  the  data  for  frame 
selection,  while  excellent  seeing  conditions  put  frame  selection  to  the  test  of  sorting 
frames  that  are  very  similar. 

3.4  Bispectrum  Parameters 

3.4-1  Phase  Spectrum  Reconstruction.  Section  2.5  discusses  the  process 
of  obtaining  the  bispectrum  estimate,  Bu{fi,f2),  from  the  detected  image  and  cor¬ 
recting  for  noise  effects.  The  sample-based  mean  bispectrum  is  an  accumulation 
of  data  from  each  of  the  frames  in  the  selection.  The  phase  spectrum  of  the  ob¬ 
ject  can  be  derived  from  the  bispectrum.  There  are  two  methods  available  in  the 
computer  simulation:  the  recursive  method  and  the  iterative  weighted  least  squares 
method.  The  recursive  method  for  reconstruction  is  the  easiest  and  most  widely 
used  [1,  12,  14,  15,  19,  24].  This  method  uses  the  fact  that  the  object  phase  spec¬ 
trum  at  {fi  -\-  f-i)  can  be  expressed  as  a  unit  amplitude  phaser,  [24] 

exp{j0o(/i  +  A)}  =  exp{i(/)o(/i)}  -I-  exp{j(j)o{f2)}  -  expO'<^B(A,  A)},  (3.2) 

where  (f)B{fi,f2)  is  the  phase  of  the  mean  bispectrum.  The  three  starting  points  for 
the  recursive  process  are  set  equal  to  zero: 

</>o(/=0)  =  0 

<;6„(±A/,0)  =  0 

</>o(0,±A/)  =  0  (3.3) 

While  <f>o{f  =  0)  is  known  to  be  zero,  zero  is  an  arbitrary  choice  for  the  other  two 
terms  in  Eq.  (3.3).  The  consequence  of  choosing  an  arbitrary  zero,  versus  a  known 
value  in  Eq.  (3.3)  is  a  hnear  shift  of  the  image  in  space  [24],  also  known  as  tilt 
error.  It  is  shown  by  Eq.  (2.17)  that  the  bispectrum  is  insensitive  to  tilt  error.  Once 
the  estimated  phase,  4>o{f))  is  reconstructed,  it  is  then  combined  with  the  estimated 
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modulus  of  the  object,  Q(/),  to  reconstruct  the  image  of  the  object.  As  a  measure  of 
performance,  the  phase  error  and  comparative  phase  error  (discussed  in  section  2.5.4) 
are  calculated  and  plotted. 

3.5  Image  Reconstruction 

The  final  image  is  composed  using  the  pseudo- Wiener  filter  described  in  section 
2.4.1  in  conjunction  with  the  modulation  transfer  function  (MTF)  filter  discussed 
in  this  section.  The  MTF  is  the  modulus  of  the  OTF  [8,  11].  The  a/SNRQ{f) 
term  in  Eq.  (refeq:wiener)  reduces  the  effects  of  high  spatial  frequency  noise  in  the 
reconstructed  image  [24].  The  results  in  this  thesis  were  generated  by  using  an  a 
value  of  0.01.  The  fc  used  in  Eq.  (3.4)  for  the  MTF  filter  in  the  computer  simulation 
is  same  as  the  calculated  OTF  radius  of  30  pixels. 

Fig.  3.2  shows  the  three  reference  star  images  used  in  the  deconvolution  pro¬ 
cess.  In  each  case,  a  bright  star  was  simulated  regardless  of  the  visual  magnitude  of 
the  satellite,  although  the  three  different  images  correspond  to  different  seeing  con¬ 
ditions.  The  reference  star  for  the  poor  seeing  conditions,  tq  =  7  cm,  is  noticeably 
more  blurred  around  the  edges  of  the  image. 

Once  the  modulus  is  obtained  from  the  pseudo-Wiener  filter  calculation,  and 
the  phase  has  been  obtained  from  the  bispectrum  calculations,  the  next  step  is  image 
reconstruction.  Before  the  final  step  of  inverse  Fourier  transforming  the  spectrum 
back  into  the  image,  another  filter  is  applied  to  correct  for  additional  high  frequency 
noise  effects.  There  are  two  filters  to  choose  from,  the  pill  box,  and  the  MTF,  or 
conical  filter.  The  MTF  filter  was  chosen  for  this  thesis  effort,  because  it  is  less  noisy 
and  tapers  off  the  high  frequency  noise  better  than  the  pill  box  filter.  The  filter  is 
applied  to  the  estimated  modulus  of  the  object  spectrum,  0(/),  obtained  from  the 
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Figure  3.2  Reference  star  images  used  in  the  pseudo- Wiener  filter  at  three  different 
Fried  seeing  levels:  a)  ro  =  7  cm  b)  r,,  =  12  cm  c)  To  =  17  cm 
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square  root  of  |0(/)p,  from  Eq.  (2.12).  It  is  defined  by  [24] 

H(i)  =  1-^.  1/1  </. 

Jc 

=  0,  |/|>/„  (3.4) 

where,  /c  is  the  OTF  radius.  The  object  spectrum  is  obtained  by  multiplying  the 
modulus,  the  phase,  and  the  filter  of  Eq.  (3.4): 

0(f)  =  d(f)H{f)  exp{Uf)}-  (3.5) 

Finally,  the  image  is  obtained  by  taking  the  inverse  Fourier  transform  of  0{f): 

o(r)  =  yO(/)exp{j7r/-r}d/.  (3.6) 


3.6  Summary 

This  chapter  provided  the  methodology  for  conducting  the  computer  simulation 
of  creating  speckle  images,  applying  frame  selection,  using  the  bispectrum  calcula¬ 
tions  for  estimating  the  phase  for  image  reconstruction,  and  applying  filters  in  the 
image  reconstruction  process.  This  chapter  also  identified  specific  parameters  used 
in  conducting  this  research.  The  final  quality  of  the  composite  image  is  determined 
by  the  power  spectrum  SNR  and  the  bispectrum  phase  error.  These  metrics  are  the 
tools  used  in  the  analysis  of  two  experiments:  the  point  source  and  the  extended 
source.  The  results  of  this  analysis  axe  presented  in  the  next  chapter. 
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IV.  Results 


4-1  Introduction 

Chapter  three  developed  the  methodology  and  introduced  the  independent 
parameters  and  performance  metrics  for  conducting  this  investigation.  This  chapter 
consists  of  the  results  of  two  frame  selection  experiments:  point  source  and  extended 
source.  The  point  source  is  typically  a  single  bright  star.  Point  source  data  is  of 
interest  to  astronomers,  and  is  used  in  most  deconvolution  techniques  for  sharpening 
images.  Satellites  are  the  extended  objects  of  primary  interest  to  the  Air  Force. 
Figure  2.1  shows  the  computer  rendering  of  a  satellite  used  as  the  extended  object 
in  this  study.  The  results  are  in  the  form  of  power  spectrum  SNR,  bispectrum  phase 
error  analysis.  The  reconstructed  images  are  also  presented. 

4.2  System  Parameters 

In  both  experiments  two  independent  parameters  are  varied:  the  Fried  seeing 
parameter  and  the  visual  magnitude.  The  Fried  seeing  parameters  are  defined  in 
section  2.7.2.  The  three  r,,  values  tested  in  this  thesis  are:  7  cm;  12  cm;  and  17  cm. 
Visual  magnitude,  m^,  is  defined  in  2.7.1.  This  investigation  explores  objects  with 
mi,  values  of  +1,  +3,  +5,  and  +7. 

4.3  Quality  Metric  Performance 

This  section  presents  a  comparison  of  the  performance  of  the  five  sharpness 
metrics  tested  in  both  the  point  source  and  extended  source  experiments.  Results  are 
given  in  terms  of  power  spectrum  SNR  and  bispectrum  phase  error.  For  illustration 
purposes,  the  vertical  scale  for  the  point  source  phase  error  diagrams  is  expanded 
much  more  than  the  extended  source. 
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4-. 3.1  Point  Source.  Figure  4.1  shows  the  radially  averaged  power  spectrum 
SNR  and  the  Fourier  phase  error  of  the  composite  point  source  simulated  with  an  Vo 
of  7  cm  and  a  of  1.  These  are  the  parameters  for  a  bright  target  in  poor  seeing 
conditions.  This  figure  plots  the  results  for  all  of  the  sharpness  metrics  discussed 
in  section  2.6.3.  The  frame  selection  rate,  FSR,  is  50%.  Figures  A.13-A.24  in  the 
appendix  are  SNR  plots,  and  figures  A.1-A.12  are  phase  error  plots  of  point  source 
images  simulated  with  other  parameters. 

When  there  are  very  poor  seeing  conditions,  the  S7  metric  proves  superior. 
This  is  because  the  Sr  metric  measures  maximum  entropy.  Recall  from  section  2.6.3 
the  Sj  metric  is 

S7  —  -  J  J ln[i{f)]i{f)  df  (4.1) 

As  seeing  conditions  worsen,  the  increased  randomness  results  in  increased  errors 
in  the  image  plane,  resulting  in  poor  quality  images.  Bright  sources  in  poor  seeing 
conditions  are  the  conditions  the  S7  metric  seems  to  work  best  in.  The  increased 
variance  due  to  the  poor  seeing  conditions  combined  with  a  high  contrast  due  to 
the  bright  source  is  a  favorable  condition  for  the  Sj  metric  to  differentiate  between 
frames.  However,  improving  the  seeing  conditions  by  simply  changing  the  value 
from  7  to  12  cm  results  in  a  dramatically  different  performance  as  can  be  observed 
in  Fig.  4.2.  In  fact,  the  results  of  the  S7  metric  were  every  inconsistent  beyond  the 
Vo  =  7cm  point.  This  inconsistency  can  be  observed  by  comparing  figures  A.3-A.6. 
Therefore,  to  produce  reliable  results,  a  more  consistent  metric  is  desired. 

The  metrics  that  performed  the  most  consistently  were  Si  and  Fig.  4.2  is 
typical  of  the  relative  performance  of  the  Ssi  and  the  Si  metrics.  Here,  they  result 
in  the  highest  SNR  and  the  lowest  phase  error.  However,  considering  the  vertical 
scale  in  this  example,  and  in  most  of  the  cases  tested,  all  of  the  metrics  performed 
about  the  same.  The  loss  in  SNR  and  the  resulting  error  achieved  by  frame  selection 
for  a  point  source,  even  in  the  worst  case,  is  insignificant.  The  worst  phase  error  is 
less  than  0.15  rad^. 
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Contrary  to  expectations,  in  most  of  the  cases  tested,  the  fidelity  defect  metric, 
Ss,  generally  resulted  in  the  lowest  SNR.  Additionally,  Ss  resulted  in  increased  phase 
error  over  the  other  metrics  in  many  of  the  scenarios.  See  also,  figures  A.2,  A. 5, 
and  A. 6.  Since  the  fidelity  defect  also  requires  a  known  source,  it  is  only  practical 
as  a  basis  for  comparison  in  an  adaptive  optics  simulation. 

4.3.2  Extended  source.  Fig.  4.3  is  representative  of  a  bright  extended 
source,  ruv  =  1,  under  average  seeing  conditions,  =  12  cm.  While  there  was 
some  loss  in  the  SNR  due  to  frame  selection,  the  metrics  all  performed  about  the 
same  for  an  extended  source.  There  was  a  slight  increase  in  the  phase  error  in 
the  mid-spatial  frequency  regions.  The  data  tended  to  fluctuate  as  the  frequency 
increased.  Figures  B.1-B.12  in  the  appendix  show  that  this  result  was  consistent  for 
all  parameters. 

4.3.3  Conclusions.  The  S'7  metric  performed  the  best  in  poor  seeing  con¬ 
ditions  for  a  bright  point  source,  but  performed  the  worst  under  improved  seeing 
conditions.  The  and  the  Ssi  metrics  performed  most  consistently,  and  performed 
well.  There  was  very  little  difference  in  metric  performance  for  an  extended  source. 
Since  Ssi  performed  consistently  well,  and  has  been  proven  superior  in  previous 
works  [5,  25],  it  will  be  used  to  represent  the  general  performance  of  frame  selection 
through  out  the  remainder  of  this  chapter. 

4.4  FYame  Selection  Performance 

4.4-^  Point  Source.  Fig.  4.4  shows  the  results  for  a  bright  point  source, 
—  1,  under  average  seeing  conditions,  ro  =  12  cm.  Here,  the  only  metric  used  is 
Ssu  but  the  frame  selection  rates  are  vaxied  from  100%  to  50%.  Fig  4.4  a)  shows 
a  gradual  decline  in  the  SNR  performance  of  frame  selection  for  a  point  source  as 
the  FSR  drops  from  100%  to  50%.  The  total  drop  in  SNR  is  about  20%,  overall. 
Figure  2.4,  is  a  graph  of  the  same  100%  data  as  in  figure  4.4,  except  the  data  is  rep- 
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(b) 

Figure  4.3  Quality  metric  performance  for  an  extended  source:  a)  SNR  and  b) 
Phase  error.  Both  cases  indicate  frame  selection  strategies  for  the  50% 
FSR.  Visual  parameters  are  ro=12  cm,  m^—l 


(a) 


Figure  4.4  Frame  selection  performance  for  a  point  source:  a)  SNR  and  b)  Phase 
error.  Both  cases  indicate  frame  selection  strategies  using  the  Ssi  per¬ 
formance  metric.  Visual  parameters  are  ro—l2  cm,  m,,=l 
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resents  the  theoretical  reduction  in  SNR  due  to  randomly  selected  data,  as  opposed 
to  using  any  sharpness  metric.  Comparing  figure  4.4  a)  to  figure  2.4  shows  that 
reducing  data  using  frame  selection  results  in  improved  SNR  over  random  selection. 

Figure  4.4  b)  shows  a  slight  improvement  in  phase  error.  According  to  the 
vertical  scale  in  b),  the  improvement  in  phase  error  at  the  mid-spatial  frequency 
region  is  less  than  one  quarter  of  a  graduation,  0.005  rad^.  Typically,  however,  the 
50%  FSR  case  is  only  as  good  or  slightly  worse  than  the  other  FSRs.  This  trend 
can  be  observed  by  comparing  the  Ssi  curves  in  figures  A.2-A.3,  in  the  appendices. 
These  figures  are  all  results  of  poor  seeing  conditions,  r-o  =  7  cm.  These  figures 
indicate  that  the  60%  FSR  performs  best.  Figure  A.4  indicates  that  the  80%  FSR 
performs  best.  In  the  succeeding  figures,  it  becomes  diflScult  to  compare  relative 
performance,  since  for  improved  seeing  conditions,  there  is  fittle  variance  in  FSR 
performance.  The  results  identified  here  are  consistent  when  other  parameters  are 
used.  Other  parameters  can  be  observed  by  referring  to  figures  A.1-A.12.  Each 
figure  has  four  plots:  one  for  eaich  FSR.  The  different  FSRs  can  be  compared  by 
selecting  a  specific  metric  curve  from  each  plot  and  observing  the  trend  as  the  FSR 
changes. 

4-4-^  Extended  source.  The  extended  source  SNR  results  are  consistent 
with  the  point  source  SNR  results.  Fig.  4.5  a)  indicates  about  a  20%  drop  in  SNR 
for  the  50%  FSR  case.  Fig.  4.5  b)  shows  that  the  phase  error  performance  of  all 
FSRs  tested  are  about  the  same,  and  the  data  tended  to  fluctuate  in  the  high  spatial 
frequency  regions.  This  result  is  consistent  with  those  using  other  parameters,  and 
can  be  observed  by  referring  to  figures  B.13-B.24. 

4.4-3  Conclusions.  Figures  4.4  and  4.5  indicate  that  there  is  a  gradual 
worsening  in  the  SNR  performance  of  frame  selection,  for  both  point  and  extended 
sources,  as  the  FSR  decreases.  This  is  consistent  with  the  statistical  nature  of 
reducing  data,  discussed  in  section  2.6.2.  Although,  the  SNR  achieved  using  the  Ssi 
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(a) 


Figure  4.5  Frame  selection  performance  for  an  extended  source:  a)  SNR  and  b) 
Phase  error.  Both  cases  indicate  frame  selection  strategies  using  the 
Ssi  performance  metric.  Visual  parameters  are  ro=12  cm,  mj,—! 
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metric  is  improved  over  that  of  randomly  selected  data.  The  phase  error  performance 
did  not  change  significantly  as  a  function  of  FSR  in  either  the  point  or  extended 
source  experiments.  In  general,  frame  selection  performance  decreased  as  the  FSR 
was  decreased.  However,  this  degradation  was  slight  enough  that  the  50%  FSR 
was  sufficient  to  illustrate  the  overall  frame  selection  performance.  For  this  reason, 
the  50%  FSR  will  be  used  throughout  the  remainder  of  this  chapter.  The  results 
obtained  by  using  the  other  FSRs  are  included  in  the  appendices. 

4.5  Performance  Under  Various  Seeing  Conditions 

The  seeing  conditions  tested  are 

•  7  cm  (poor  seeing  conditions) 

•  12  cm  (average  seeing  conditions) 

•  17  cm  (excellent  seeing  conditions). 

SNR  gain  and  comparative  phase  error  difference  plots  are  used  in  this  section 
to  show  relative  frame  selection  performance  between  cases  where  seeing  conditions 
are  varied.  The  SNR  gain  is  discussed  in  section  2.6.2,  and  phase  error  difference  is 
discussed  in  section  2.5.4.  Each  data  set  in  the  graphs  represents  relative  performance 
based  on  the  50%  FSR. 

4-5.1  Point  Source.  To  create  figure  4.6,  a  bright  point  source  in  varying 
seeing  conditions  is  simulated.  Here,  the  FSR  is  fixed  to  50%,  and  the  sharpness 
metric  used  is  Ssi-  Each  of  the  three  ro  values  are  compared.  The  complete  set 
of  SNR  gain  plots  are  in  figures  A.25-A.36.  These  figures  include  aU  the  quality 
metrics  and  all  of  the  FSRs.  To  determine  the  relative  performance  of  different 
seeing  conditions,  three  figures  of  different  r,,  values,  but  the  same  values  must 
be  compared.  For  example,  figure  4.6  a)  displays  the  SNR  gain  data  from  the  Ssi 
curve  in  the  a)  plots  from  figures  A. 25,  A.29,  and  A. 33.  The  SNR  gain  in  Fig.  4.6  a) 
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shows  a  slight  decline  in  frame  selection  performance  as  seeing  conditions  get  worse. 
In  nearly  every  case  studied,  the  SNR  gain  curves  showed  a  slight  improvement  only 
in  regions  of  low  frequency,  and  a  slight  decline  elsewhere. 

The  complete  set  of  comparative  phase  error  plots  are  in  figures  A.37-A.48. 
The  comparative  phase  error  is  introduced  in  section  2.5.4.  Figure  4.6  b)  dis¬ 
plays  the  comparative  phase  error  data  from  the  Ssi  curve  in  the  a)  plots  from 
figures  A.37,  A.41,  and  A.45.  The  comparative  phase  error  plot  in  b)  also  shows 
that  performance  declines  when  the  seeing  conditions  drop  below  =  12  cm.  Al¬ 
though  the  plot  in  b)  has  such  a  small  scale  that  there  is  really  very  little  significant 
difference  in  performance.  The  lowest  point  in  the  curve  is  about  0.1  rad^. 

4-5.2  Extended  Source.  In  figure  4.7,  a  bright  extended  source  in  varying 
seeing  conditions  is  simulated.  Here,  the  FSR  is  fixed  to  50%,  and  the  sharpness 
metric  used  is  Ss\.  Each  of  the  three  To  values  are  compared.  The  complete  set  of 
SNR  gain  plots  are  in  figures  B.25-B.36.  Figure  4.7  a)  displays  the  SNR  gain  data 
from  the  Ssi  curve  in  the  a)  plots  from  figures  B.25,  B.29,  and  B.33. 

The  SNR  gain  plot  in  Fig.  4.7  a)  indicates  a  20%  drop  in  performance  for 
the  50%  frame  selected  case.  It  also  shows  that  the  performance  is  about  the  same 
for  all  seeing  conditions.  The  complete  set  of  comparative  phase  error  plots  are  in 
figures  B.37-B.48. 

Figure  4.7  b)  displays  the  comparative  phase  error  data  from  the  Ss\  curve 
in  the  a)  plots  from  figures  B.37,  B.41,  and  B.45.  The  phase  error  differences  in 
Fig.  4.7  b)  show  there  is  no  significant  change  in  performance  whether  or  not  frame 
selection  is  used,  at  any  of  the  values  tested.  In  nearly  every  case  studied,  the 
low  spatial  frequency  region  of  the  spectrum  has  almost  zero  phase  error  difference 
between  the  frame  selected  case  and  the  non-frame  selected  case.  As  the  frequency 
increases,  the  phase  error  difference  oscillates  about  the  zero  axis. 
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Comparative  Phase  Error  Difference  (rad'^2)  Comparative  SNR  Gain 


(b) 

Figure  4.6  Seeing  conditions  for  a  point  source:  a)  SNR  gain  curves  and  b)  Phase 
error  difference  curves.  Both  cases  indicate  frame  selection  strategies 
for  the  50%  FSR,  using  the  Ssi  performance  metric.  m^=l 
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Comparative  Phase  Error  Difference  {rad'^2)  Comparative  SNR  Gain 


rO  =  1 7  cm 


Figure  4.7  Seeing  conditions  for  an  extended  source:  a)  SNR  gain  curves  and  b) 
Phase  error  difference  curves.  Both  cases  indicate  frame  selection  strate¬ 
gies  for  the  50%  FSR,  using  the  performance  metric.  771^=1 
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4.5.3  Conclusions.  The  Vo  values  tested  only  affect  frame  selection  perfor¬ 
mance  for  point  sources.  As  seeing  conditions  improve,  so  does  performance.  In  the 
case  of  extended  sources,  the  effect  of  changing  To  is  insignificant.  These  results  are 
consistent  with  the  results  using  the  other  m,,  parameters  tested. 

4.6  Performance  Under  Various  Target  Brightness  Levels 

The  brightness  levels  tested  are  m„  values  of  -f-1,  -1-3,  -(-5,  and  +7. 

4.6.1  Point  Source.  In  figure  4.8,  a  varying  brightness  point  source  in 
average  seeing  conditions  is  simulated.  To  determine  the  relative  performance  of 
varying  brightness  levels,  four  figures  from  the  appendix  containing  the  same  Tq 
value  but  different  values  may  be  compared.  Figure  4.8  a)  displays  the  SNR  gain 
data  from  the  Ssi  curve  in  the  a)  plots  from  figures  A.29,  A.30,  A.31,  and  A. 32. 
Fig.  4.8  a)  shows  there  is  a  very  slight  drop  in  frame  selection  SNR  performance  as 
the  target  brightness  level  decreases.  The  worst  case  tested  was  =  -t-7,  which 
had  about  a  20%  reduction  in  SNR  performance. 

Figure  4.8  b)  displays  the  comparative  phase  error  difference  data  from  the 
Ssi  curve  in  the  a)  plots  from  figures  A.41,  A.42,  A.43,  and  A.44.  Figure  4.8  b)  also 
indicates  a  drop  in  phase  error  performance  as  target  brightness  decreases. 

4.6.2  Extended  Source.  In  figure  4.9,  a  varying  brightness  extended  source 
in  average  seeing  conditions  is  simulated.  Figure  4.9  a)  displays  the  SNR  gain  data 
from  the  Ssi  curve  in  the  a)  plots  from  figures  B.29,  B.30,  B.31,  and  B.32.  Fig.  4.9 
shows  there  is  a  very  slight  drop  in  frame  selection  SNR  performance  for  aU  target 
brightness  levels,  especially  at  the  low  frequencies.  There  is  no  significant  reduction 
in  performance  in  the  phase  error  for  any  of  the  brightness  levels.  Once  again, 
performance  begins  to  oscillate  as  the  frequency  increases,  however,  staying  very 
close  to  the  performance  levels  of  the  non-frame  selected  case. 
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Comparative  Phase  Error  Difference  (rad*2)  Comparative  SNR  Gain 


Spatial  Frequency  (normalized) 


(b) 

Figure  4.8  Brightness  levels  for  a  point  source:  a)  SNR  gain  curves  and  b)  Phase 
error  difference  curves.  Both  cases  indicate  frame  selection  strategies 
for  the  50%  FSR,  using  the  Ssi  performance  metric.  ro=12  cm 
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Comparative  Phase  Error  Difference  (rad''2)  Comparative  SNR  Gain 


(b) 


Figure  4.9  Brightness  levels  for  an  extended  source:  a)  SNR  gain  curves  and  b) 
Phase  error  difference  curves.  Both  cases  indicate  frame  selection  strate¬ 
gies  for  the  50%  FSR,  using  the  Ssi  performance  metric.  ro=12  cm 


4-6.3  Conclusions.  Lower  visual  magnitude  values  tend  to  reduce  perfor¬ 
mance  for  point  sources.  The  visual  magnitude  values  tested  have  no  significant 
effect  on  the  performance  of  frame  selection  for  the  extended  sources.  These  results 
are  consistent  with  other  To  parameters  tested. 

4.7  Image  Reconstruction  Results 

Figure  4.10  shows  the  reconstructed  image  of  an  extended  source  a)  before 
frame  selection  and  b)  after  a  50%  FSR  frame  selection  was  imposed.  There  is  no 
noticeable  difference  in  final  image  quality  whether  100%  or  50%  of  the  frames  are 
retained  for  post  processing. 

4.7.1  Conclusions.  A  50%  FSR  frame  selection,  using  the  551  metric  does 
not  reduce  the  visual  quality  of  speckle  images. 

4.8  Summary 

Frame  selection  is  an  effective  way  to  reduce  data  by  as  much  as  50%,  thus 
reducing  the  expense  and  time  of  post-processing  speckle  image  data,  by  nearly 
a  factor  of  two.  The  quality  sharpness  metrics  varied  little  in  performance.  In 
situations  where  the  seeing  conditions  are  very  poor,  rg  <  12  cm,  the  maximum 
differential  sharpness  metric  {S^)  is  superior.  The  spectral  maximization  metric 
{Ss\)  was  chosen  to  represent  frame  selection  performance  because  of  its  consistency 
to  perform  well  over  a  wide  variety  of  conditions.  Although,  frame  selection  SNR 
performance  was  reduced  by  about  20%  for  point  sources,  the  phase  error  did  not 
increase  significantly.  The  benefit  gained  by  a  50%  data  reduction  outweighs  the 
20%  drop  in  performance.  Furthermore,  when  it  comes  to  extended  sources,  the 
primary  Air  Force  interest,  frame  selection,  using  the  metric  with  a  50%  FSR 
performed  as  well  as  the  100%  case  over  all  ro  and  values  tested. 
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(b) 

Figure  4.10  Typical  extended  source,  reconstructed  using  a)  100%  of  the  frames  b) 
50%  frame  selection,  using  the  spectral  maximization  [Ssi )  metric,  via 
a  1  meter  telescope  at  500  kilometers,  with  seeing  conditions  of  ro=07 
cm  and  m„=l. 
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V.  Conclusions  and  Recommendations 


This  investigation  began  due  to  the  need  to  devise  an  effective  post-processing 
technique  that  does  not  rely  on  expensive  adaptive  optics.  Frame  selection  needed 
to  be  extended  to  uncompensated  speckle  imaging.  A  complete  understanding  of  the 
capability  of  this  technique  is  critical  to  its  application  to  U.  S.  Air  Force  imaging 
problems.  This  thesis  determines  whether  or  not  frame  selection  processing  enhances 
the  speed  and  efficiency  of  data  management  for  image  processing,  without  significant 
loss  in  quality.  This  quality  is  measured  in  terms  of  power  spectrum  SNR  and 
bispectrum  phase  error.  This  chapter  presents  a  summary  of  the  accomplishments 
documented  in  this  thesis  and  provides  recommendations  for  further  research  in  this 
area. 

5.1  Conclusions 

This  thesis  makes  the  following  conclusions: 

1.  The  five  quality  sharpness  metrics  all  performed  about  the  same  for  both  point 
and  extended  sources. 

2.  Frame  selection  using  the  Ssi  metric  performed  well  in  all  circumstances  for 
both  point  and  extended  sources.  The  50%  FSR  resulted  in  a  20%  reduction 
in  SNR  performance  for  the  point  source,  but  no  significant  change  in  phase 
error.  For  extended  sources,  there  was  no  significant  change  in  SNR  or  phase 
error  performance. 

3.  Increasing  the  Fried  seeing  parameter,  Tg,  improves  image  quality,  but  has 
little  effect  on  frame  selection  performance.  Decreasing  Vg  had  the  effect  of 
decreasing  frame  selection  SNR  performance  for  point  sources,  but  no  effect 
on  phase  error.  Decreasing  Vg  had  no  effect  on  frame  selection  performance  for 
extended  sources. 
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4.  Increased  target  brightness,  in  the  form  of  decreased  visual  magnitude, 
improves  image  quality,  but  has  little  effect  on  frame  selection  performance. 
Increasing  m„  had  the  effect  of  decreasing  frame  selection  SNR  performance 
for  point  sources,  but  no  effect  on  phase  error.  Increasing  rrii,  had  no  effect  on 
frame  selection  performance  for  extended  sources. 

5.  A  50%  frame  selection  using  the  Ssi  sharpness  metric,  for  the  extended  source, 
results  in  a  composite  image  that  has  the  same  visual  quality  as  the  100%  case. 
Satellites  are  the  extended  sources  are  of  primary  interest  to  the  U.  S.  Air  Force. 
The  benefit  gained  from  frame  selection  is  a  factor  of  two  reduction  in  data 
required  to  produce  the  composite  satellite  image.  A  50%  reduction  in  data 
translates  to  nearly  double  the  processing  speed. 

5.2  Recommendations  for  Further  Research 

1.  Frame  selection  should  be  fully  demonstrated  on  real  data  collected  with  ade¬ 
quately  short  exposure  times.  This  data  could  be  obtained  from  an  Air  Force 
facility  such  as  the  Starfire  Optical  Range  (SOR)  3.5  meter  telescope,  or  the 
Air  Force  Maui  Optical  Site  (AMOS)  1.6  meter  telescope. 

2.  The  effects  of  exposure  time  on  image  quality  should  be  fully  demonstrated  on 
speckle  image  processing.  The  simulation  package  used  in  this  thesis  does  not 
take  into  account  the  effect  of  increasing  the  integration  or  exposure  times,  r,  of 
each  image  data  frame.  In  an  operating  environment,  there  exists  an  optimum 
value  for  r.  This  optical  parameter  has  both  positive  and  negative  effects.  This 
thesis  has  only  taken  into  account  the  positive  effect  of  increasing  r,  that  is, 
increasing  the  total  number  of  photons  per  image.  This  by  itself  translates  into 
increased  target  brightness,  which  is  a  measure  of  image  quality.  However,  this 
single  assumption  neglects  the  fact  that  increasing  r  has  the  negative  effect  of 
increasing  the  blurring  caused  by  turbulence  and  target  motion,  thus  decreasing 
the  image  quality. 
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3.  Other  metrics,  involving  the  unbiased  image,  or  the  power  and  phase  spectrums 
should  be  explored.  Since  each  metric  studied  so  far  only  examines  a  single 
characteristic  of  image  features,  a  new  metric  might  contain  aspects  of  more 
than  one  of  the  metrics  studied  in  this  thesis. 

4.  Since  speckle  imaging  has  not  been  shown  to  improve  as  a  result  of  frame  selec¬ 
tion,  closer  attention  needs  to  be  paid  to  how  far  the  performance  is  reduced.  A 
comparison  of  frame  selected  data  to  randomly  selected  data  deserves  further 
consideration. 
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Appendix  A.  Simulation  Results  for  Point  Source 

The  appendices  include  the  complete  collection  of  data  results  for  all  the  quality 
sharpness  metrics  discussed  in  section  3.3.2.  FSR’s  of  100,  80,  70,  60,  and  50  %  are 
all  included.  These  are  the  results  of  two  experiments:  the  point  source,  and  the 
extended  source.  In  each  experiment,  two  visual  parameters  were  varied:  seeing 
conditions,  Vo,  and  visual  magnitude,  m^,.  The  results  of  these  experiments  are 
presented  in  four  parts:  1)  the  bispectrum  phase  error;  2)  the  power  spectrum  SNR; 
3)  the  power  spectrum  SNR  gain  (with  vs.  without  frame  selection);  and  4)  the  4) 
comparative  phase  error  plots. 

This  appendix  contains  the  results  of  the  point  source  experiment. 
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A.l  Point  Source  Bispectrum  Phase  Error 

This  section  contains  the  phase  error  results  for  the  point  source  experiment. 
On  each  page,  the  four  FSR  strategies  are  presented,  containing  the  results  of  the 
five  metrics  tested.  A  given  set  of  visual  parameters,  ro  and  m„,  are  identified  below. 
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Radially  Averaged  Fourier  Phase  Errc 
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A. 2  Point  Source  Power  Spectrum  SNR 

This  section  contains  the  power  spectrum  SNR  results  for  the  point  source 
experiment.  On  each  page,  the  four  FSR  strategies  are  presented,  containing  the 
results  of  the  five  metrics  tested.  A  given  set  of  visual  parameters,  and  mj,,  are 
identified  below. 
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Figure  A.  17  Power  spectrum  SNR, 
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A. 3  Point  Source  Power  Spectrum  Gain  Plots 

This  section  contains  the  power  spectrum  SNR  gain  plots  for  the  extended 
source  experiment.  Gain  is  defined  as  the  ratio  of  the  power  spectrum  achieved  by 
using  frame  selection  to  that  without  frame  selection.  On  each  page,  the  four  FSR 
strategies  axe  presented,  containing  the  results  of  the  five  metrics  tested.  A  given 
set  of  visual  parameters,  r^  and  are  identified  below. 
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Figure  A. 32  Power  spectrum  gain,  point  source,  =  7,  To  =  12  cm.  FSR  =  a) 
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Figure  A.36  Power  spectrum  gain,  point  source,  =  7,  ro  =  17  cm.  FSR  =  a) 
50%  b)  60%  c)  70%  d)  80% 
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A. 4  Point  Source  Comparative  Phase  Error  Difference  Plots 

This  section  contains  the  phase  error  difference  plots  for  the  point  source  exper¬ 
iment.  Phase  error  difference  is  defined  as  the  difference  in  the  phase  errors  between 
the  cases  using  frame  selection  to  that  without  frame  selection. 

The  relative  phase  error  performance  of  a  case  with  frame  selection  is  compared 
to  that  without  by 

<^a(/*)  =  100%  ~  FS^  (■^•1) 

where  <f)e{f)ioo%  represents  the  phase  error  of  the  100%  case,  and  4>e{f)FS  represents  a 
frame  selected  case.  Areas  in  the  spectrum  where  performance  is  improved  as  a  result 
of  frame  selection,  will  have  positive  values.  Conversely,  where  the  performance 
declines  as  a  result  of  frame  selection,  (pA  values  are  negative. 
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Figure  A.42  Comparative  phase  error  difference,  point  source,  m,,  =  3,  Tq  =  12  cm. 
FSR  =  a)  50%  b)  60%  c)  70%  d)  80% 
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Figure  A.47  Comparative  phase  error  difference,  point  source,  rriv  =  5,  Tq  =  17  cm. 
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Appendix  B.  Simulation  Results  for  Extended  Source 


This  section  contains  the  results  of  the  extended  source  experiment.  The  results 
of  this  experiment  is  presented  in  four  parts:  1)  the  bispectrum  phase  error;  2)  the 
power  spectrum  SNR;  3)  the  power  spectrum  SNR  gain  (with  vs.  without  frame 
selection);  and  the  compaiative  phase  error  plots. 
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B.l  Extended  Source  Bispectrum  Phase  Error 

This  section  contains  the  phase  error  results  for  the  extended  source  experi¬ 
ment.  On  each  page,  the  four  FSR  strategies  are  presented,  containing  the  results 
of  the  five  metrics  tested.  A  given  set  of  visual  parameters,  rg  and  are  identified 
below. 


B-2 


Averaged  Fourier  Phase 


(c)  (d) 

Figure  B.2  Bispectrum  phase  error,  extended  source,  —  3,  ro  =  7  cm.  FSR 
a)  50%  b)  60%  c)  70%  d)  80% 
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Figure  B.3  Bispectrum  phase  error,  extended  source,  rUi,  —  5,  To  —  7  cm.  FSR 
a)  50%  b)  60%  c)  70%  d)  80% 
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B.2  Simulation  Results  for  Extended  Source 

This  section  contains  the  power  spectrum  SNR  results  for  the  extended  source 
experiment.  On  each  page,  the  four  FSR  strategies  are  presented,  containing  the 
results  of  the  five  metrics  tested.  A  given  set  of  visual  parameters,  Xo  and  are 
identified  below. 
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Figure  B.20  Power  spectrum  SNR,  ext' 
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Figure  B.23  Power  spectrum  SNR,  extended  source,  =  5,  Vo  =  17  cm.  FSR 
a)  50%  b)  60%  c)  70%  d)  80% 


B.3  Extended  Source  Power  Spectrum  Gain  Plots 

This  section  contains  the  power  spectrum  SNR  gain  plots  for  the  point  source 
experiment.  Gain  is  defined  as  the  ratio  of  the  power  spectrum  achieved  by  using 
frame  selection  to  that  without  frame  selection.  On  each  page,  the  four  FSR  strate¬ 
gies  are  presented,  containing  the  results  of  the  five  metrics  tested.  A  given  set  of 
visual  parameters,  ro  and  m^,  are  identified  below. 
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B.4  Extended  Source  Comparative  Phase  Error  Difference  Plots 

This  section  contains  the  phase  error  diiference  plots  for  the  extended  source 
experiment.  Phase  error  difference  is  defined  as  the  difference  in  the  phase  errors 
between  the  cases  using  frame  selection  to  that  without  frame  selection. 

The  relative  phase  error  performance  of  a  case  with  frame  selection  is  compared 
to  that  without  by 

^a(/*)  =  (pe{f)^lQ0%  ~  PS')  (B-1) 

where  <^e(/)ioo%  represents  the  phase  error  of  the  100%  case,  and  4>e{f)FS  represents  a 
frame  selected  case.  Areas  in  the  spectrum  where  performance  is  improved  as  a  result 
of  frame  selection,  will  have  positive  values.  Conversely,  where  the  performance 
declines  as  a  result  of  frame  selection,  (pA  values  are  negative. 

On  each  page,  the  four  FSR  strategies  are  presented,  containing  the  results  of 
the  five  metrics  tested.  A  given  set  of  visual  parameters,  ro  and  m^,  are  identified 
below. 
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